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Abstract 

We provide a model-independent framework for the analysis of the radiative 
.B-meson decays B — > A*7 and B — > pry. In particular, we give a systematic 
discussion of the various contributions to these exclusive processes based on the 
heavy-quark limit of QCD. We propose a novel factorization formula for the con- 
sistent treatment of B — > V'y matrix elements involving charm (or up-quark) 
loops, which contribute at leading power in ^-qcd / m B to the decay amplitude. 
Annihilation topologies are shown to be power suppressed. In some cases they are 
nevertheless calculable. The approach is similar to the framework of QCD fac- 
torization that has recently been formulated for two-body non-leptonic B decays. 
These results allow us, for the first time, to compute exclusive b — > s(d) r y decays 
systematically beyond the leading logarithmic approximation. We present results 
for these decays complete to next-to-leading order in QCD and to leading order 
in the heavy-quark limit. Phenomenological implications for various observables 
of interest are discussed, including direct CP violation, and isospin and U-spin 
breaking effects. 
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1 Introduction 



The radiative transitions b — > 37, b — > d'y are among the most valuable probes 
of flavour physics. Proceeding at rates of order G"j?a, they are systematically 
enhanced over other loop-induced, non-radiative rare decays, which are propor- 
tional to G 2 F a 2 . In fact, the Cabibbo-favoured b — > S7 modes belong to the small 
number of rare decays that are experimentally accessible already at present. The 
inclusive branching fraction has been measured to be 

B(B -> X s7 ) = (2.96 ±0.35) ■ 10" 4 (1) 

combining the results of fl], ^|, ^ . The branching ratios for the exclusive channels 
have been determined by CLEO and more recently also by BABAR || and 
BELLE 101: 



B(B° -> K*°j) = \ (4.39 ±0.41 ±0.27) • 10- 5 [§ (2) 



(4.55 ±0.70 ±0.34) ■ 
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10~ 5 


i 


(3.76 ±0.86 ±0.28) 


• 10~ 5 
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(3.89 ±0.93 ±0.41) 


• 10" 5 
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B(B+^K* + -i) I (3.89 ±0.93 ±0.41) • 10~ 5 § ^ 

On the theoretical side, the flavour- changing neutral current (FCNC) reactions 
b — > s(d) , ~f are characterized by their high sensitivity to New Physics and by the 
particularly large impact of short- distance QCD corrections || 0, || II- Con- 
siderable efforts have therefore been devoted to achieve a full calculation of the 
inclusive decay b —>■ s 7 at next-to-leading order (NLO) in renormalization group 
(RG) improved perturbation theory [II], [XT], [T2|] (see |nj for recent reviews). 

Whereas the inclusive mode can be computed perturbatively, using the fact 
that the 6-quark mass is large and employing the heavy-quark expansion, the 
treatment of the exclusive channel B — > K*^ is in general more complicated. In 
this case bound state effects are essential and need to be described by nonper- 
turbative hadronic quantities (form factors). The basic mechanisms at next-to- 
leading order were already discussed previously for the B — > amplitudes |14 . 



However, hadronic models were used to evaluate the various contributions, which 
did not allow a clear separation of short- and long-distance dynamics and a clean 
distinction of model-dependent and model-independent features. 

In this paper we present a systematic analysis of the exclusive radiative decays 
B — > (V = K*, p) in QCD, based on the heavy quark limit nib ^qcd- We 
shall establish factorization formulas for the evaluation of the relevant hadronic 
matrix elements of local operators in the weak Hamiltonian. Factorization holds 
in QCD to leading power in the heavy quark limit. This result relies on arguments 
similar to those used previously to demonstrate QCD factorization for hadronic 
two-body modes of the type B — > nn [T5|, |T5| . 



This framework will allow us to separate perturbatively calculable contribu- 
tions from the nonperturbative form factors and universal meson light-cone dis- 
tribution amplitudes (LCDA) in a systematic way. This includes the treatment of 
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loop effects from light quarks, in particular up and charm. Such loop effects are 
straightforwardly included for the inclusive decays b — > s(d)^f. For the exclusive 
modes, however, the effects from virtual charm and up quarks have so far been 
considered to be uncalculable "long-distance" contributions and have never been 
treated in a model independent fashion. 

Finally, power counting in KqcD/ m b implies a hierarchy among the possible 
mechanisms for B — > V'j transitions. This allows us to identify leading and 
subleading contributions. For example, weak annihilation contributes only at 
sub leading power in the heavy quark limit. 

Within this approach, higher order QCD corrections can be consistently taken 
into account. We give the B — > decay amplitudes at next-to- leading order 
(NLO). After including NLO corrections the largest uncertainties still come from 
the B —* V form factors, which are at present known only with limited precision 
(~ ±15%), mostly from QCD sum rule calculations fll7| . The situation should 
improve in the future with the help of both lattice QCD (TB| and analytical 
methods based on the heavy-quark and large-energy limits |21| . 



Despite the more complicated theoretical situation of the exclusive channels in 
comparison to the inclusive decays, the goal of obtaining a better understanding 
of the exclusive modes is very well motivated. This is because the exclusive decays 
are easier to investigate experimentally, especially in the difficult environment of 
hadron machines as the Fermilab Tevatron or the LHC at CERN. In any case the 
systematic uncertainties, both experimental and theoretical, are very different for 
B — > V^i and b — > s(d)7- A careful study of the exclusive modes can therefore 
yield valuable complementary information in testing the Standard Model. 

The remainder of this paper is organized as follows. In section |2| we collect 
basic expressions and present the structure of the QCD factorization formulas 
for B — > V7 matrix elements in general terms. Sections |3| and |] give a detailed 
discussion of B — > K*^ and B — > fry, respectively. Phenomenological applications 
are studied in section ^. Finally, section |6] contains our conclusions. 

2 Basic Formulas 

The effective Hamiltonian for b — > 57 transitions reads 



U ff = — V A (s) 

V " p=u,c 



C\Q\ + C2Q2 + X/ CiQi 



i=3,...,8 



(4) 



where 



A? } = V; s V pb (5) 



The operators are given by 



Qi = {sp)v-A{pb)v-A (6) 

Ql = {SiPj) V ~A{pjbi)v-A (7) 
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Qz = 


(sb) v ^ A J2(Q<l)v-A 
q 


Q 4 = 


q 


Q 5 = 


(s6) y _A^(gg)y+A 

q 


Q 6 = 


q 


Qi = 


g 

m b Sicr^(l + 75)6^ 


Qs = 


^m 6 5^(1 + 75)^ ^ 



(9) 
(10) 

(11) 

(12) 
(13) 

The most important operators are Q\ 2 an d the magnetic penguin operator Q-j. 
Note that in our notation the numbering of Q\ 2 is reversed with respect to the 
convention of p2j, i.e. Cf ] {M w ) = 1 and Cf\M w ) = 0. The sign conventions 



for the electromagnetic and strong couplings correspond to the covariant deriva- 
tive = + ieQfA^ + igT a A a ^. With these definitions the coefficients Cr,a 
are negative in the Standard Model, which is the choice generally adopted in the 
literature. The effective Hamiltonian for b — > d^ is obtained from (j^TB]) by the 
replacement s — > d. The Wilson coefficients Cj in (|j) are known at next-to-leading 
order 



The most difficult step in computing the decay amplitudes is the evaluation of 
the hadronic matrix elements of the operators in @). A systematic treatment can 
be given in the heavy-quark limit. We will argue that in this case the following 
factorization formula is valid 



(V-y(e)\Qi\B) 



Jo 



(14) 



where e is the photon polarization 4-vector. Here F B ^ V is a B — > V transition 
form factor, and $£, are leading twist light-cone distribution amplitudes of 
the B meson and the vector meson V, respectively. These quantities are universal, 
nonperturbative objects. They describe the long-distance dynamics of the matrix 
elements, which is factorized from the perturbative, short-distance interactions 
expressed in the hard-scattering kernels T/ and T? 1 . The QCD factorization 
formula ( ITS ) holds up to corrections of relative order KqcD/ m b- 

For Qi the factorization formula (|14}) is trivial. The matrix element is simply 
expressed in terms of the standard form factor, T/ is a purely kinematical function 
and the spectator term T*, 1 is absent. An illustration is given in Fig. |]. The 
matrix element reads 

& r 

(V(k,r])j(q, e)\Q 7 \B) = -—m b c v F v [e^ xp e^ v k x q p + i(e ■ rjk ■ q - e ■ krj ■ q) 

(15) 

where Cy = 1 for V = K*, p and Cy = l/y/2 for V = p°. The B — > V form 
factor Fy is evaluated at momentum transfer q 2 = 0. Our phase conventions 
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Figure V. Contribution of the magnetic penguin operator Q-j described by B — > V 
form factors. All possible gluon exchanges between the quark lines are included 
in the form factors and have not been drawn explicitly. 





Figure 2: 0{ot s ) contribution to the hard-scattering kernels T/ from four-quark 
operators Qi. The crosses indicate the places where the emitted photon can be 
attached. 



coincide with those of [O, 23]. In particular we have Fy > 0, e = —1, and 



the phases of V (with flavour content qq') and B are such that 

(V(k,7])\qa IJil/ q , \0) = -i(%k„ - tyk^fv (16) 
(0\u W b\B(p)) = +if BP , (17) 

with positive fg, fy- In the leading logarithmic approximation (LO) and to 
leading power in the heavy-quark limit, Q-j gives the only contribution to the 
amplitude of B — > Vj. 

The matrix elements of the four-quark operators Qi (and of Qq) start con- 
tributing at 0(a s ). In this case the factorization formula becomes nontrivial. The 
diagrams for the hard-scattering kernels T/ are shown in Fig. ||] for Qi, ■ ■ ■ ,Qe 
and in Fig. |^ for Q s . The non- vanishing contributions to T? 1 are shown in Fig. 

The diagrams in Figs. - ^ represent the complete set of contributions from 
<5x,2 an d Qs to leading order in the heavy-quark limit and at 0(a s ). The first term 
in the factorization formula ([TJ]) is given by the form factor from ( |T5| ) and the 
hard-scattering kernels T i . The latter are calculable functions of m q /mf,, where 
m q is the mass of the internal quarks in the loop diagram. For the second term 
in (0), where the spectator quark is involved, we need light-cone distribution 



4 



X r >< 



Figure 3: C?(a s ) contribution to the hard-scattering kernels T/ from chromomag 
netic penguin operator Q 8 . 




b 



^< 5, d 



Figure 4: 0(a s ) contribution to the hard-scattering kernels Tf 1 from four-quark 
operators Qi (left) and from Q 8 . 



amplitudes (wave functions) for B mesons and light vector mesons. In the case 
of the B meson we have to leading power Jl6|, |20 

(0\b(0)u(z)\B(p)) = l Jf^ + mb ) l5 [ 1 d£e- i ^*-[<S> B1 (0+ (18) 

4 Jo 

with n = (1,0,0,-1), chosen to be parallel to the 4-momentum of the vector 
meson. The functions $1,2(0 describe the distribution of light-cone momentum 
fraction £ = l+/p+ of the spectator quark with momentum I inside the B meson. 
Here light-cone components of four-vectors v are defined by 



v± 



V° ± V 3 

~7T 



(19) 



The wave functions are highly asymmetric with £ = 0(AQCD/ m b)- They are 
normalized as 

[ 1 d^ Bl (o = i /V<M£) = o (20) 

Jo Jo 
The first negative moment of $bi(£), which will be needed below, can be para- 
metrized by a quantity \ B = O(Aqcd), i-e. 



1 d ^ m ^ = — 
ii £ X B 



(21) 



In B — > Vj decays the vector meson is transversely polarized. The leading- 
twist and leading-power distribution amplitude for light vector mesons with trans- 
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Figure 5: Annihilation contribution to B — > V'j decay. The dominant mechanism 
is the radiation of the photon from the light quark in the B meson, as shown. 
This amplitude is suppressed by one power of Ag OT /m b , but it is still calculable 
in QCD factorization. Radiation of the photon from the remaining three quark 
lines is suppressed by (Aqcd/^) 2 for operators <2i,2- 



verse polarization, is defined by 



(V(k, V )\q'(z)q(0)\0} = ^cx^K f dv e^Q^v) (22) 

4 Jo 

Here and in the following we use the short-hand notation 

v = 1 - v (23) 

for light-cone variables. 

The light-cone wave function $j_ has an expansion in terms of Gegenbauer 
polynomials C^/ 2 (2v — 1) 

$ ± (y) = Q V {1 



n=l 



(24) 



where C/ (x) = 3x, C 2 (x) = |(5x 2 — 1), etc. The Gegenbauer moments 
a^(/z) are multiplicatively renormalized. They vanish logarithmically as the scale 
H — > oo. In this limit $j_ reduces to its asymptotic form $j_(f) = 6vv, which 
often is a reasonable first approximation. The remaining leading-twist light-cone 
wave functions for light vector mesons, $H , g { l ] and fl7|, H, do not contribute 
at leading power if the mesons are transversely polarized. 

The form factor Fy, /b/^b, and the light-cone wave function /y$±(f) are 
the nonperturbative quantities required to describe B — > at next-to-leading 
order in QCD. ($B2(0 does not contribute.) 

There are further mechanisms that can in principle contribute to B — > 
decays. One possibility is weak annihilation, depicted in Fig. |5|. In this case the 
leading-power projection onto the meson V in (^) vanishes, because the trace 
over an odd number of Dirac matrices is zero. A non- vanishing result arises from 
the projector 

(V(k, V )\q lu q'\0) = f v m vVu (25) 
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Figure 6: Other contributions that are power-suppressed in the heavy-quark limit. 



which, however, is suppressed by one power of compared to (|IED, (E2p 

for transverse polarization rj v (/y, fy, my ~ Aq C d, k ~ m^). This results in a 
corresponding power suppression of weak annihilation. The dominant contribu- 
tion comes from the diagram shown in Fig. |5]. Here the photon is emitted from 
the light-quark constituent in the B meson, which leads to a quark propagator 
scaling as 1/Aqcd- This is in contrast to the remaining three possible diagrams 
with a quark propagator ~ 1/rrib. The latter contributions are therefore even 
stronger suppressed, ~ (Aqcd/"^) 2 relative to the leading B — > amplitude. 

Despite its power suppression, the dominant annihilation amplitude can be 
computed within QCD factorization. This is because the colour-transparency 
argument applies to the emitted, highly energetic vector meson in the heavy- 
quark limit |nj. A similar observation was already made in PS|. On the other 
hand, we disagree with the claim made in this paper that the 0(a s ) correction 
to Fig. H] would vanish identically in the chiral limit and to leading-twist order. 
This claim is based on the observation that the projector in (p2|) gives zero when 
applied to a current with an odd number of Dirac matrices. However, by the 
same argument, the diagram in Fig. [| would vanish even at leading order in 
a s , which is not the case. The proper treatment of the 0(a s ) correction should 
employ the subleading-power projections related to the wave functions $y, g±\ 

g^f 1 Jl7|, 124}! ) corresponding to the use of (p5|) at 0{oP s ). This correction has not 
yet been computed, but we see no a-priori reason to expect a vanishing result. 

Since weak annihilation is a power correction, we will content ourselves with 
the lowest order result (O(a s )) for our estimates below. In particular, we shall 
include the annihilation effects from operators to estimate isospin-breaking 
corrections in B — > p7 decays. The reason for including this class of power 
corrections is that they come with a numerical enhancement from the large Wilson 
coefficients Ci j2 (Ci ~ 3|C 7 |) and are not CKM suppressed. Instead, a CKM 
suppression of annihilation effects occurs for B —>■ K*^ and these contributions 
are thus very small in this case. 

Finally, Fig. |6| displays the remaining diagrams with a power suppression. 
The penguin-annihilation diagrams (left) are power suppressed in a way similar 
to the ordinary annihilation topologies. The spectator diagrams from Qg (right) 
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lead to amplitudes that are either manifestly power-suppressed (photon emitted 
from the quark line in the upper right), or that are superficially of leading power, 
but vanish when the leading-order projections are performed (photon emitted 
from either of the quarks forming the B meson). Recall, however, that photon 
emission from the light-quark line from Q 8 is a leading-power effect (see Fig. |]). 



3 B — > IP7 

In the case of B — > K*^ the component of the Hamiltonian (^) proportional to 
\ u is strongly CKM suppressed (|A U /A C | ~ 0.02) and has only a minor impact on 
the decay rate. It is essentially negligible, but will be included later on for com- 
pleteness. Throughout this work we shall neglect the contribution from the QCD 
penguin operators Q3, . . . , Qq, which enter at 0{a s ) and are further suppressed 
by very small Wilson coefficients. We note that to 0(a s ) the matrix element of 
Q 2 is zero because of its colour structure. The amplitude for B — > K*^ then 
reads «&) = (K^Q^Bj) 

A(B - K*j) = ^|A« (C 7 (Q 7 ) + C^Ql) + C S (Q 8 )) (26) 



The leading-order matrix element (Qr) is given in (15) with V = K*. At sub 



leading order in a s the matrix elements (Qi.s) need to be computed from the 
diagrams in Figs. |2| - [|. The result for the diagrams in Figs. |2] and |||, which enter 
the hard-scattering kernels T{, T/, can be infered from [IT|. In these papers the 



diagrams were computed to obtain the matrix elements for the inclusive mode 
b — > S7 at next-to-leading order. In this context Figs. ^ and |3] represented the 
virtual corrections to the inclusive matrix elements of Qi and Q$. In our case 
they determine the kernels T( and T/. As required for the consistency of the 
factorization formula these corrections must be dominated by hard scales ~ m;, 
and hence must be infrared finite. This is indeed the case. Re-interpreted as the 
perturbative hard-scattering kernels for the exclusive process, the results from 

im p!y 

(Qi,8>' = (Q 7 )^G 1 , 8 (27) 
where Cp = (N 2 — 1)/(2N), with N = 3 the number of colours, and 

G 1 (s c ) = -^\n-^ + 9l (s c ) (28) 
27 m b 

G 8 = ^ln-^ + g 8 (29) 
3 m h 

833 20m 8tt 2 



9i[s) = ~i62-^r + — s 

2 

+ 9 



3/2 

48 + mix - 5vr 2 - 2ivr 3 - 36C(3) + (36 + 6m - 9tt 2 ) In s 



8 



+ 9 

1 

+ 



+ (3 + Gin) ln 2 s + ln 3 s 



18 + 2vr 2 - 2i?r 3 + (l2 - 6tt 2 ) In s + 6ivr ln 2 s + ln 3 s 
9 + 112Z7T - 14vr 2 + (182 - A8in) In s - 126 ln 2 s 



27 

11 27T 2 2/77 

"3" ~ ~9 



+ 



Here we denote 



s 3 (30) 
(31) 



(32) 



We now turn to the mechanism where the spectator participates in the hard 
scattering. 

To find the correction for (Qi) we compute the first diagram in Fig. |], using 
the light-cone projectors in (|T8|) and (p2|) . We obtain 



(Qi) u = (Qr) 



An 



H x (s c ) 



with 



Hi{s) 



3NF v m 2 B Jo ( Jo 

The hard-scattering function h(u,s) is given by 



dv h(v, s)$j_(f ) 



/i(w, s) 



4s 



/ 'I — — ; : 

L 2 is the dilogarithmic function 



+ L 2 



1 + 



U — l.S + tg 



2 
u 



(33) 
(34) 

(35) 



L,(x) = - dt 



* ln(l - t) 



(36) 



The function /i(w, s) is real for w < 4s and develops an imaginary part for u > As. 
At small values of u it has the expansion 



h(u, s) 

It is also regular for s — ► 



1 1 ^ / « 2 

— + — M + — 

6s 45s" 5 V 



h{u,0) 



u 



(37) 



(3f 



The function h(v, s c ) is displayed in Fig. |7[ 

The correction to (Q&) from the hard spectator interaction comes from the 
second diagram in Fig. |j. One finds 



(Qs) n = (Q 7 )^^H 8 

An 



(39) 




Figure 7: The hard-scattering kernel h(v, s c ) as a function of v. 



where 



4VT 2 hft f' ^BliQ /- 1 $lQQ 



(40) 



The gluons in Fig. f| transfer a momentum of order ~ ^AQCDf^b- There- 
fore we set a s = a s (fih) in (5D and fl39|) . For our numerical analysis we shall use 
Hh = y/Xhjl with Aft = 0.5 GeV and \i = 0(m&). 

Finally, we can combine these results and write, including also the up-quark 
contribution, 



A(B -> iT 7 ) 



Gp 
V2 



E^i s) «?(^7) 



(X*7lQr|5) 



(41) 



where, at NLO 



a?(y 7 )=C 7 + aa( y F (gi(/i)Gi(ap)+C 8 MG 8 ) 



+ 



47T 

An 



{CMH^ + CMH^) 



(42) 



Here the NLO expression for CV has to be used, while the leading order values 
are sufficient for C\ and C§. The explicit formulas for the Wilson coefficients can 



be found in 12 . 
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The scale dependence of the matrix element of Qj is reflected in the running 
of the product of 6-quark mass and form factor which is explicitely given as 

(m b ■ F v ) \p] = {m b ■ F v ) [m b ] (l - In (43) 

This dependence on the scale p has to be taken into account when the residual 
scale dependence of physical quantities is investigated. 

From the amplitude in (|4T|) the branching ratio is obtained as 



B{B^K* 1 ) = r B G> ^ 1 (l-^) \Y.^^1)?\F^ (44) 



4 B 



PI 



For the decay B — > pj both sectors of the effective Hamiltonian have the same 
order of magnitude and have to be included. It is straightforward to translate the 
expressions from the previous section to this case. The amplitude can be written 

as 



A{B -» 

where a^iV^) is given in 



\ Gf 



L V 



{pi\Qi\b) 



(45) 



The branching fraction becomes 



B (B -> pry) = r B 



G 2 F am B ml 
32^ 



fl-^) 3 |EAf ^(pi)\ 2 cl\F p 

\ 111 B J p 



(46) 



The rate for the CP-conjugated mode B — > p^ is obtained by replacing A^ 
\w*. We may then consider the CP asymmetry 



A C p{pi) 



T{B -> pry) - T(B -> pry) 
r(B^p 1 ) + F(B^p 1 ) 



(47) 



A non-vanishing CP asymmetry appears at 0(a s ) only. Expanding Acp in 
a s and using the improved Wolfenstein parametrization for the CKM elements 
we obtain 



Acp(pl) 



IA 



(<9|2 

2?7 



|C 7 | 2 



(4? 



;i - p) 2 + f7 2 4tt C< 



' X S Im (Gi(s c ) + Hi(s c ) - d(0) - #i(0)) 



As we have discussed above, weak annihilation from the leading operators Qi,2 
contributes to the B — > p7 amplitude only at C^Ag^/m?,), but is enhanced by 
large Wilson coefficients. These effects can be calculated in QCD factorization. 
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To lowest order in a s we obtain from the diagram in Fig. |5|, using the projection 
inflU) 

A ann (B-^p-j) = ( ^L\fa 1 b u { p - 1 \Q 7 \B-) (49) 
A ann (B° -> p° 7 ) = ^Ai^MAlQrlS ) (50) 

ai,2 = Ci, 2 + — C 2 ,i (51) 



Here 



(evaluated in leading logarithmic approximation) and 



2ir 2 QufBf P m P b _ lj 
F p m\\ B d 2 



6u = ^ 1 P &<* = ^ (52) 



Recalling that F p ~ m b 3 ^ 2 , fs ~ 1//2 in the heavy-quark limit, we note that 
b u ~ KQcol^b- This shows explicitly the power suppression of weak annihilation. 
The ratio of bd to b u is —Qd/Qu — 1/% where the minus sign comes from the 
relative sign between the up-quark and down-quark components of the p° wave 
function (only the up-quarks produce the p° in the annihilation process, while 
only the down-quarks are relevant in (Q7)). 

The annihilation components are included in the decay amplitudes by substi- 
tuting 

a" — > aj+ b u ai for B~ — > p~7 (53) 
-> al} + b d a 2 for B° -> p° 7 (54) 

The annihilation contribution for 5 _ — > K*~ , -f is similar to the one for B~ — > 
p~7, with obvious replacements p — > K* in 6 U . Its impact is very small for B — > 
K*~^ because of the strong CKM suppression of the up-quark sector. It will be 
included in our numerical analysis for defmiteness. Within our approximations 
weak annihilation gives no contribution to B° — > K* 0r y. 

We conclude this section with a few general remarks on the role of power 
corrections in B — > Vj decays. Let us first repeat that the annihilation effect 
from operator Qi gives a numerically important power correction, because it 
comes with a relative enhancement factor of |Ci/CV| ~ 3. This leads to a 30% 
correction in the amplitude of the charged mode B~ — > p~7- Note that the 
parameter b u in (|o"2|), describing the generic effect of the annihilation term, is only 
about 10%, which is consistent with a hQCD/ m b correction of canonical sizeQ. In 

^^After the preprint version of the present article had appeared, further isospin-breaking 
power corrections were investigated in ^6|. A sizeable effect of 11% from penguin annihilation 
related to Qq was identified. This is larger than naively expected, but still substantially smaller 
than the effect of b u (for the p channel) , which is presumably the largest power correction. The 
Qe annihilation contribution is still calculable, while other terms of the same order are much 
smaller numerically. 
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contrast to nonleptonic modes with pseudoscalar mesons in the final state, no 
chirally enhanced power corrections from the light-cone expansion of meson wave 
functions arise in our case. Finally, power corrections can also come from the 
loop effects with up- and charm quarks, whose leading-power contributions were 
computed in (|33|) . These power corrections correspond to the region of integration 
where the gluon becomes soft, that is v = 0(AQ C D/ m b)- Their contribution 
is nonperturbative and cannot be calculated in the hard-scattering formalism. 
Nevertheless, the expression in ( ]3"B| ) can be used to read off the scaling behaviour 
of these power corrections in the heavy-quark limit. For the charm loop the 
kernel approaches a constant ~ ml/m 2 c in the endpoint region as shown in ([37l). 
Taking into account the linear endpoint suppression of the wave function 
the integral in (|33| ) over the region v ~ hQCD/ m b therefore contributes a term 
of order {Aq C D/ m b) 2 x (jnb/ m c) 2 = (Aqcd/^c) 2 - It is interesting that we thus 
recover the power behaviour of soft contributions in the charm sector first pointed 
out in ^7|. This was discussed for the inclusive decay b — > 57 in J27|, and 



for the exclusive mode B — > K*j in [29]. Numerically this correction is very 



small (~ 3% in the decay rate). A similar consideration applies to the up-quark 
sector. In this case the endpoint behaviour of the kernel is singular fl3"8|), which 
now leads to a linear power suppression of the form KQcoj^b- This coincides 
with the scaling behaviour derived in [28 in the context of the inclusive process. 

In view of these remarks, we see no indication for large uncalculable power 
corrections. Further studies on this issue are desireable, but would go beyond the 
scope of the present analysis. 



5 Phenomenology 

In this section we present numerical results for various observables of interest 
for the phenomenology of B — > V'-f decays, based on the NLO QCD expressions 
derived in this paper. Our choice of input parameters is summarized in Table |l|. 
Our default choice for the CKM angle 7 is 58 deg. 

We begin with the numerical result for the NLO QCD coefficient a^(K*'~f) as 
a typical example. For central values of all parameters, at fi = nib, and displaying 
separately the size of the various correction terms, we find 

a c 7 (K*j) = -0.3221 +0.0113 -0.0820 - 0.0147i -0.0144 -0.0109* 
Cf° ACf LO T/ 8 -contribution ^-contribution 

= -0.4072 - 0.02561 (55) 

We note a sizable enhancement of the leading order value, dominated by the T 1 - 
type correction. This feature was already observed in the context of the inclusive 
|TTj| . A complex phase is generated at NLO, where the T^-corrections and 



case m 



the hard-spectator interactions (T 11 ) yield comparable effects. 

The net enhancement of a-j at NLO leads to a corresponding enhancement 
of the branching ratios, for fixed value of the form factor. This is illustrated 
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Table 1: Summary of input parameters. 



CKM parameters and coupling constants 


V us 


Vcb 


\v ub /v ct 


| 


MS 


a 


Gp 


0.22 


0.041 


0.085 ±0.025 


225 MeV 


1/137 


1.166 x 10- 5 GeV" 2 


Parameters related to the B mesons 


m B 


h 




T B + 


T B 


5.28 GeV 


180 MeV 


(350 ± 150) MeV 


1.65 ps 


1.56 ps 


Parameters related to the K* meson 


|1 


Zl 






m K * 


<* 


af 


f K * 


0.38 ±0.06 


185 MeV 


894 MeV 


0.2 


0.04 


230 MeV 


Parameters related to the p meson 
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F P 


J P 


m p 


a{ 




f P 


0.29 ±0.04 


160 MeV 


770 MeV 





0.2 


200 MeV 


Quark and W-boson masses 


m b (m h ) 


m c {m b ) 


m t,pole 


M w 


4.2 GeV 


(1.3 ±0.2) GeV 


174 GeV 


80.4 GeV 



Table 2: Predictions for branching ratios and CP asymmetries with the errors 



from the individual input uncertainties. 





B(B° ^ir*° 7 )[10- 5 ] 


B(B- ^p- 7 )[10- 6 ] 




central 


7.09 


1.58 


9.89 


Fr* 


+2.32/-1.99 






F P 




±0.42/-0.37 


+0.12/-0.13 




+0.67/-1.00 


+0.07/-0.19 


±5.25/-2.40 


Xb 


+0.38/-0.15 


±0.26/-0.09 


±0.35/-0.28 


m c 


±0.38/-0.43 


+0.12/-0.12 


+1.40/-1.52 


\v ub /v cb \ 


±0.05/-0.05 


+0.15/-0.13 


±3.48/-3.30 
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H [GeV] ^ [GeV] 

Figure 8: Dependence of the branching fractions B(B° — > K* 0/ ~f) and B(B~ 
p~7) on the renormalization scale p at leading and next-to-leading order. 



in Fig. where we show the residual scale dependence for B[B — > A'*°7) and 
B[B~ — > p~7) at leading and next-to-leading order. 

The sensitivity of the B — > ^~*°7 and — > p~7 branching ratios, and of 
the CP asymmetry Acp{pi) to variations in the relevant input parameters are 
summarized in Table 2. The uncertainty of the branching fractions is currently 
dominated by the form factors Fk*, F p . The values for B(B~ — > K*~^f) are very 
close to those for B(B — > K* 0r y) with the main shift coming from the different 
lifetimes. For central input parameters we find B(B~ — > K*~j) = 7.45 ■ 10~ 5 
compared to B(B -> K*°~/) = 7.09 ■ 10" 5 . The CP asymmetry A C p{K*^) is 
typically -0.5%. 

Taking the sizable uncertainties into account, the results for B — > K*^ in 
Table 2 are compatible with the experimental measurements in (0) and (|3|), even 
though the central theoretical values appear to be somewhat high. 

The direct CP asymmetry, which is substantial for the pj modes, is much less 
dependent on the form factors. Here the largest theoretical uncertainty comes 
from the scale dependence. This is to be expected because the direct CP asym- 
metry is proportional to the perturbative strong phase difference, which arises 
at 0(a s ). Unknown power corrections could have some impact on the predic- 
tion. B — > pj also depends sensitively on fundamental CKM parameters, such 
as |T4&/Kb| and 7, and can thus in principle serve to constrain the latter quan- 
tities once measurements become available. This is further illustrated in Figs. 
and [U| where the dependence on 7 is shown for B(B~ — > p~j) and Acp(p^), 
respectively. 
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Figure 9: The branching fraction B(B 
7 at leading and next-to-leading order. 



p 7) as a function of the CKM angle 



Further interesting observables are the isospin breaking quantities 

A - T{B+ - p+l) 1 (56) 

A - T{B ~ - P ' 7) 1 (57) 

A(P7) = (58) 

Within our approximations, isospin breaking is generated by weak annihilation. 
Isospin breaking was already discussed in pOj , partially including NLO correc- 
tions. Our results, including the power-suppressed annihilation effects at leading 
order in QCD, together with the complete NLO expressions at leading power, 
are displayed in Fig. [TTj. We remark that our sign of A(p7) differs from the one 
found in []30] . 



Another application of our results concerns an estimate of U-spin breaking 
effects in B — > Vj decays. U-spin symmetry, the symmetry of strong interactions 
under exchange of d and s quarks, has been advocated as a tool to control hadronic 
uncertainties in tests of the Standard Model pi] , [J^, Defining 



AB(B - 
AB(B 



-> pi) 



B{B + 
B(B + 



K*+i) - B(B- 
p+ 7 )-B(B- 



- JT- 7 ) 
P'l) 



(59) 
(60) 



16 



16 



14 



12 



i 10 



a. 

U 



|j.=m h /2 




20 40 60 80 100 120 140 160 180 

y[deg] 

Figure 10: The CP asymmetry Acp{pi) as a function of the CKM angle 7 for 
three values of the renormalization scale \x = m&/2, m& and 2mi,. 



the quantity 



AB(B -> iT T ) + A£(£ -> p 7 ) = 



(61) 



in the limit of U-spin symmetry and within the Standard Model. This has been 
discussed in PT[ and was considered in more detail in [K2] . Using our expressions 



and central values for all parameters we find 



AB(B - 
AB(B 



iT 7 ) 
-> P7) 



-7- 10" 7 
4- 1(T 7 



(62) 
(63) 



where we have chosen the CKM angle 7 = tt/2, which maximises the effects. 
The two quantities indeed have opposite signs, but their sum only partly cancels, 
leaving a U-spin breaking remainder of —3 • 10~ 7 . This effect is almost entirely 
due to the difference (Fk* — F p ). For form-factor values different from those in 
Table 1, the U-spin breaking effect would approximately be rescaled proportional 
to (Fk* — F p ). For our choice the sum of the two asymmetries in (|6i|) is of the 
same order of magnitude as the individual asymmetries. This example quantifies 



the limitations of the relation (|BTD as a Standard Model test. 
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20 40 60 80 100 120 140 160 180 

y[deg] 



Figure 11: The isospin-breaking asymmetry A(fry) as a function of the CKM 
angle 7 at leading and next-to-leading order. 

6 Conclusions 

In this paper we have discussed a systematic and model-independent framework 
for the exclusive radiative decays B — > V'-f based on the heavy-quark limit. This 
allowed us to compute the decay amplitudes for these modes consistently at next- 
to-leading order in QCD. 

An important conceptual aspect of this analysis is the interpretation of loop 
contributions with charm and up quarks, which come from leading operators in 
the effective weak Hamiltonian. We have argued that these effects are calculable 
in terms of perturbative hard-scattering functions and universal meson light-cone 
distribution amplitudes. They are 0(a s ) corrections, but are leading power con- 
tributions in the framework of QCD factorization. This picture is in contrast to 
the common notion that considers charm and up-quark loop effects as generic, un- 
calculable long-distance contributions. Non-factorizable long-distance corrections 
may still exist, but they are power-suppressed. 

Another important feature of the NLO calculation are the strong interaction 
phases, which are calculable at leading power. They play a crucial role for CP 
violating observables. 

We have seen that weak-annihilation amplitudes are power-suppressed, but 
can be numerically important for B — > because they enter with large coef- 
ficients. These effects also turn out to be calculable and were included in our 
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phenomenological discussion at leading order in QCD. 

Finally, we have presented a numerical analysis of important observables, 
including the branching ratios of B — > K* r y 1 B — > /ry, and CP asymmetries and 
isospin-breaking in B — > pj decays (see sec. 0). Currently, the large uncertainties 
in the B — > V form factors are still an important limitation, but the situation can 
be systematically improved. In particular, our approach allows for a consistent 
perturbative matching of the nonperturbative form factor to the short-distance 
part of the amplitude. Our formalism can also be applied to other radiative rare 
B decays, such as B s — > Vj, V = 0, K* , or Bd — > ivy (see e.g. p5f for a discussion 
of these modes). 

The improved theoretical understanding of B — > Vy decay observables streng- 
thens the motivation for still more detailed experimental investigations, which will 
contribute significantly to our knowledge of the flavour sector. 
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Note Added 



After completion of this work we received the manuscript J55 ] discussing B —>■ fq 
decays at next-to-leading order in QCD. The calculation uses the framework 



of |15|, [U| and of p0| , and is similar in spirit to the analysis presented here. 
We disagree, however, with several of the main results given in fl3~5|| . First, the 
expressions corresponding to our functions i?i(0) and H 8 are missing a factor of 



—4 and —1, respectively, in f5H) |. Furthermore, the charm contribution Hi(s c ) is 
approximated by i?i(0). As can be seen from our results, this approximation is 
not suitable. It gives a real part with the wrong sign and misses a complex phase, 
which is important for Acpipj) (changing the asymmetry in Table 2 from 9.9% 
to 13.2%). Finally, we find Hi(s c ) to be free of infrared divergences, also for non- 
vanishing charm-quark mass, contrary to P5| . This is important conceptually, 



because such a divergence could not be reabsorbed into the form factor and would 
signify a breakdown of the factorization approach. 
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